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ABSTRACT 


The largest geomagnetic storms of solar cycle 24 so far occnrred on 2015 
March 17 and Jnne 22 with Dgt minima of —223 and —195 nT, respectively. 
Both of the geomagnetic storms show a mnlti-step development. We examine 
the plasma and magnetic held characteristics of the driving coronal mass ejec¬ 
tions (CMEs) in connection with the development of the geomagnetic storms. A 
particnlar effort is to reconstrnct the in sitn structnre nsing a Grad-Shafranov 
techniqne and compare the reconstruction results with solar observations, which 
gives a larger spatial perspective of the source conditions than one-dimensional 
in situ measurements. Key results are obtained concerning how the plasma and 
magnetic held characteristics of GMEs control the geomagnetic storm intensity 
and variability: (1) a sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta 
scenario are proposed for the multi-step development of the 2015 March 17 and 
June 22 geomagnetic storms, respectively; (2) two contrasting cases of how the 
GME hux-rope characteristics generate intense geomagnetic storms are found, 
which indicates that a southward hux-rope orientation is not a necessity for a 
strong geomagnetic storm; and (3) the unexpected 2015 March 17 intense geo¬ 
magnetic storm resulted from the interaction between two successive GMEs plus 
the compression by a high-speed stre am from b e hind, w hich is essentially the 
“perfect storm” scenario proposed by iLin et al.l fj2014al. i.e., a combination of 
circumstances results in an event of unusual magnitude), so the “perfect storm” 
scenario may not be as rare as the phrase implies. 
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Subject headings: shock waves — solar-terrestrial relations — solar wind — Snn: 
coronal mass ejections (CMEs) 


Introduction 


A topic of increasing interest to space weather is how the plasma and magnetic field 
characteristics of coronal mass ejections (CMEs) resnlt in geomagnetic storm activity, in 
particnlar those intense events. The sonthward magnetic field and speed of CMEs at the 
Earth have received the most attention, becanse their cross prodnct, the dawn-dnsk electric 
field, controls the rate of the solar wind energy conpling to the terrestrial magnetosphere 
fjDungevlIlQGlf) . However, it is still not clear how the ejecta speed and sonthward magnetic 
field work together to achieve a snstained, enhanced dawn-dnsk electric field and how they 
lead to the variability of geomagnetic storms. 

The sonthward magnetic field is often fonnd within the ejecta reaching the Earth in the 
form of an interplanetary CME (ICME) with a preceding shock. This nsnally leads to a clas¬ 
sic geomagnetic storm seqnence: a sndden commencement generated by the shock, a main 
decrease phase cansed by the ejecta’s sonthward magnetic field, and then a recovery phase. 
In addition to the driver gas, the sheath region betwee n the shock and ICME can also be geo¬ 


effective (e.g., iTsnrntani et al.lll988l: Ihin et al.ll2008bl) as both the sheath speed and sonth- 


invoked to explain the two-step development of geomag 
Complex eiecta resnlting from interactions between CME 

netic storms (Kamide et ah 

1998). 

3 (I 

3nrlaga et ah 2001, 20021) can be 

verv geo-effective owing to their pro 

onged dnrations (e 

s-, 

Farrngia & Berdichevskv 

2004; 

Zhang et ah 2007: Lngaz & Farrngia 

2014 

Mishra et ah 

201^. Thev could also cans 

e two- 

mplex 

step geomagnetic storms (Farrngia et al.l 2 

Ood: Ihin et ah 

2014b). A special case of co 


ejecta i s the interaction of a preceding ejecta with an overtaking shock fe.g.. lLin et al.ll2012 


2014bl: iMostl et ahl 120121: iHarrison et al.ll2012l: IWebb et al.l 120131 ). The shock enhances the 


pre-existing sonthward magne tic field inside t he ejecta, an idea fo r increased geo-effectiveness 
dating back several decades (IBnrlagalllQQll: I Vandas et al.Nl997j) . A statistical analysis in¬ 


dicates that 19 ont of 49 shocks propagating inside ICMEs are associated with an intense 
geomagnetic storm (minimnm Dgt < —100 nT nLngaz et al.ll2015l) . A recent stndy combining 
remote-sensing and in sitn observati ons snggests a “p erfect storm” scenario for the generation 
of an extreme space weather event flLin et al.ll2014al) : preconditioning of the npstream solar 
wind by an earlier CME pins in-transit interaction between later two closely lannched CMEs, 
in order to have an exceptionally high solar wind speed and nnnsnally strong ejecta magnetic 
fields at 1 AU. This, again, emphasizes the crncial importance of CME-CME interactions 
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for space weather. 


On 2015 March 17 and June 22 the Earth underwent an intense geomagnetic storm with 
the minimum D^t of —223 and —195 nT, respectively. These are the largest geomagnetic 
storms of solar cycle 24 so far. T hey occurred in the declining phase of the sol ar cycle, a phe¬ 
nomenon that is not uncommon (IGopalswamv et al.ll2005l: iKilpua et al.ll2015l ). We provide a 
timely analysis of the driving CMEs, in an attempt to identify the plasma and magnetic field 
characteristics controlling the geomagnetic storm intensity and variability . Anoth er focus of 
this Letter is to test the “perfect storm” scenario proposed bv iLiu et al.l (l2014all : whether 
it is a rare coincidence or if it happens more frequently than what the phrase suggests. We 
examine the solar wind signatures and their connections with the development of the geomag¬ 
netic storms, complemented with the r nodeling of the D^t index using two empirical form ulae 
based on the solar wind measurements ( Burton et ah 19751: O’Brien fc McPherron 200oh. W e 


also use a Grad-Shafranov (GS) technique ( Han fc Sonnerupl 19991: Hu fc Sonnerunl 2002 1. 


which has been va lidated by well separated multi-spacecraft measurements flLiu et al.ll2008al: 
Mostl et al.l 1200911 . to reconstruct the in situ IGME structure. The GS method can give a 
cross section as well as flux-rope orientation without prescribing the geometry. In conjunc¬ 
tion with solar observations it pr ovides a larger spatial perspective of IGMEs than one¬ 
dimensional in situ measurements fiLiu et al.l 1201011 . These efforts are key to understanding 
how the plasma and magnetic held characteristics of GMEs are connected with the intensity 
of geomagnetic storms as well as their variability. 


The 2015 March Event 


Tracing back to the Sun, the drivers of the 2015 March 17 geomagnetic storm were 
two interacting GMEs on March 15 (Figure 1, left). The second GME (GME2) had a 
maximum speed of about 1100 km s“^ and was associated with a long-duration G9.1 hare 
from AR 12297 (S22°W25°) that peaked at 02:13 UT on March 15. The hrst GME (GMEl) 
occurred on March 14 and had a speed of about 350 km s“^. It was likely associated with a 
G2.6 hare from the same active region (S21°W20°) that peaked around 11:55 UT on March 
14. A hrst impression from the coronagraph images is that GMEl was largely propagating 
southward while GME2 had a major component heading west. One may expect that the 
Earth would encounter the hank of the ejecta, so this would not raise the alarm for a major 
geomagnetic storm. Another fact that also makes the occurrence of an intense geomagnetic 
storm surprising is that the associated hares are relatively weak. In this sense, the occurrence 
of the 2015 March 17 intense geomagnetic storm is similar to the form ation of the 2012 July 
23 super solar storm that impacted STEREO A flLiu et ahl l2014all . Without white-light 
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observations from STEREO accurate CME kinematics cannot be obtained. 


Figure 2 shows the in situ signatures observed at Wind. A shock passed Wind at 04:01 
UT on March 17 and caused the sudden commencement of the geomagnetic storm. It is 
difficult to unambiguously connect the in situ signatures with the coronagraph images with¬ 
out wi de-angle imaging observa tions from STEREO. However, application of an empirical 
model flGopalswamv et al.ll2000[) with the CME speed of 1100 km s“^ gives a predicted ar¬ 
rival time of 23:59 UT on March 16 at Wind (0.99 AU from the Sun), which is only 4 hours 
earlier than the observed shock arrival. As can be seen from the fi gure, two ICMEs (or 
flux ropes) are identified. Our interpretation is different from those of iKataoka et ahl (120151) 
and Gopalswamv et ah ( 20151) who identify a single, shorter ICME interval from the data 
(although with different durations). The reason that we believe there are two ICMEs is as 
follows. First, there are multiple rotations in the magnetic field components whose polarities 
change twice in the shaded data intervals. Obviously these features cannot be explained by 
a single flux rope. Second, our ICME intervals are an outcome of the GS reconstruction, 
which is sensitive to the chosen boundaries. Despite the magnetic field fluctuations, both 
of the ICMEs can be reconstructed fairly well (see description below). Third, the interpre¬ 
tation of two ICMEs is consistent with what the solar observations indicate (see Figure 1 
and discussions below). One may argue for a single ICME interval (say, from 12:58 UT on 
March 17 to 03:22 UT on March 18) based on the depressed proton temperature, a signature 
often used to identify ICMEs. A reasonable GS reconstruction, however, cannot be obtained 
for the interval and its variations. Given the presence of the compression by a high-speed 
stream from behind and CME-CME interactions, the criterion of a low proton temperature 
for identifying ICMEs may not be valid in the current case. 

Within the ICME intervals the maximum magnetic field strength is about 33 nT while 
the southward component reaches —25 nT. These are not small magnetic fields considering 
only the flanks of the CMEs are encountered (see Figure 1). Interactions between these two 
ICMEs may have inhibited their expansion, a mech anism to create s trong ejecta magnetic 
fields as we have seen from the 2012 July 23 event (iLiu et al.ll2014a() . Also note the high¬ 
speed stream compressing ICME2 from behind, which may help maintain a strong ejecta 
magnetic field and a relatively high speed as well. The D^t profile indicates a two-step 
geomagnetic storm sequence with a global minimum of —223 nT. The first dip is produced 
by the southward magnetic field component in the sheath region behind the shock, while 
the second one results from the southward fields within the two ICMEs that last about 12 
hours. Given the presence of a preceding shock and two interacting ICMEs, the two-step 
development of the geomagnetic st orm can be classified a s a sh eath-ejecta-ejecta scenario. 
The modeled Dgt index using the lO’Brien fc McPherroiJ (120001) formula (minimum —170 
nT) generally agrees with actual Dgt measurements but underestimates the global minimum. 
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The iBurton et al.l fll975f) scheme gives a deeper global minimum (—280 uT) and a shallower 
recovery phase than measured. 


Further information on how the plasma and magnetic field characteristics control the 
geomagnetic storm activity is obtained from the GS reconstruction, as shown in Figure 3. 
The reconstructions give a right-handed flux rope structure for both of the ICMEs, as can 
be judged from the transverse helds along the spacecraft trajectory together with the axial 
direction. Their axis orientations are almost opposite to each other: an elevation angle 
of about 33° and azimuthal angle of about 256° (in RTN coordinates) for ICMEl, and 
an elevation angle of about —18° and azimuthal angle of about 92° for ICME2. These 
low inclinations are consistent with the slightly tilted neutral line (N. Gopalswamy, private 
communication) and hlament channel (M. Temmer, private communication) associated with 
the active region. IGMEl has a larger elevation angle (33°), which may help explain the 
encounter in spite of the largely southward propagation direction of GMEl (Figure 1). The 
impact of IGME2 may be accounted for by its lower elevation angle (—18°), although GME2 
had a major sectio n propagating westward. In addition, the angular sp an of an IGME can be 
60° or even larger (IBothmer fc Schwennlll998l: iRichardson et al.ll200211 . so it is not surprising 
that both of the IGMEs hit the Earth. It is also likely that the GME-GME interactions, 
the following high-speed stream and even the surrounding coronal magnetic held structures 


may have changed the propagation directions of both GMEs fe.g.. iGopalswamv et ah 


Znccarello et al.ll2ni2l : Ihngaz et al.ll2ni2l : [ 


Mdstl et al.lbniht iKay et al.ll2ni5l ). 


2009 


Given those relatively low inclination angles, the geomagnetic storm was mainly caused 
by the azimuthal magnetic field components of the flux ropes rather the axial components. 
A particularly interesting feature shown by the cross section of IGMEl is that the azimuthal 
field component (with a maximum value of about 30 nT) is much larger than the axial com¬ 
ponent (maximum value of only 8 uT). Therefore, if the flux rope were vertically orientated 
with the axis pointing southward the geomagnetic storm would be much weaker since the 
southward held would be less. This is contrary to the common belief that a southward point¬ 
ing hux rope favors a strong geomagnetic storm. IGME2 has comparable axial and azimuthal 
magnetic held components. Another prominent feature is the vortices visible in the cross 
section of IGME2, reminiscent of the Kelvin-Helmholtz instability. This is probably a result 
of the fast stream interacting with IGME2 from behind (see Figure 2). 


The 2015 June Event 


In 2015 June AR 12371 exhibited elevated activity, somewhat sirn ilar to AR 11429 in 
2012 March ( Liu et ah 2013, 2014c : Wang et ah 2014; Sun et ah 2015). The active region 
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produced a CME of about 1200 km s“^ associated with an M3.0 flare from N13°E45° that 
peaked at 17:36 UT on June 18, a CME of about 1300 km s“^ associated with an M2.0 flare 
from N12°E13° peaking at 01:42 UT on June 21, a CME of about 1000 km s“^ associated 
with an M6.5 flare from N13°W05° peaking at 18:23 UT on June 22, and another one of 
about 1700 km s“^ associated with an M7.9 flare from N10°W42° peaking at 08:16 UT on 
June 25. All these CMEs impacted the Earth. The June 21 CME appeared as a single halo 
event in the coronagraph images (see Figure 1, right), and a near head-on collision with the 
Earth was expected. 


The corresponding in situ signatures at Wind are displayed in Figure 4. A cluster of 
shocks passed Wind at 16:05 UT on June 21, 05:02 UT and 18:00 UT on June 22, and 13:12 
UT on June 24, respectively. The ICME boundaries are determined from the magnetic held 
in conjunction with the proton temperature and density. The hrst shock (SI) seemed driven 
by the June 18 CME, and the second one (S2) was likely associated with a CME from June 
19. No driver signatures are observed at Wind for these two shocks, presumably owing to 
the largely eastward and southward propagation directions of the June 18 and 19 CMEs 
respectively (not shown here). The ICME and its preceding shock (S3) were produc ed by 
the June 21 CME; again, application of the empirical model (IGopalswamv et al.ll2000h with 
the CME speed of 1300 km s“^ yields a predicted arrival time of 17:02 UT on June 22 at 
Wind (1.02 AU from the Sun), which is only 1 hour earlier than the observed S3 arrival. The 
fourth shock (S4) that was overtaking the ICME at 1 AU was associated with the June 22 
CME. A series of dips in the magnetic held strength are observed inside the ICME, suggestive 
of the presence of current sheets. This signature is possibly due to the heliospheric current 
sheet cutting through the ejecta, which may lead to a chain of small hux ropes within the 
ICME (see below). 

The Dst prohle shows a multi-step geomagnetic storm with a global minimum of —195 
nT. The hrst dip is produced by the huctuating southward held component upstream of S3 
(likely owing to amplihcation by the two preceding shocks), the second one by the southward 
held in the sheath downstream of S3 (further enhanced by S3), and the major dip by the 
southward held in the hrst hatched interval inside the ejecta. The southward held in the 
second hatched interval only creates a battening of the D^t value, perhaps because of the 
extremely low density. Note that the solar wind density upstream of the ejecta is signifi¬ 
cantly enhanced by the three preceding shocks. This high density may feed the plasma sheet 


of the magnetosphe re, which in turn helps intensify the ring current flFarrugia et al.l 12006 


Lavraud et al.l 120061) . Given the presence of more than one preceding shocks and a single 
IGME, the multi-step development of the geomagnetic storm can be classihed as a sheath- 
sheath-ejecta scenario. Again, the modeled Dgt index using the lO’Brien &: McPherronl fj2000[) 
formula (minimum —174 nT) underestimates the global minimum but reproduces the recov- 
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ery phase fairly well, whereas application of the iBurton et ahl fjl975[l scheme gives a larger 
global minimum (—241 nT) than measured. 


Figure 5 shows the reconstructed cross sections of two small flux ropes identihed inside 
the ICME. Here we call them small flux ropes rather than ICMEs, in order to distinguish 
from the 2015 March 17 case. Specifically, the June 22 event is a single ejecta instead of 
multiple ICMEs. These small flux ropes may have formed from the interaction between the 
CME and the heliospheric current sheet. The reconstructions yield a left-handed structure 
for both of the flux ropes and similar axis orientations: an elevation angle of about —77° 
and azimuthal angle of about 236° for the hrst one (FRl), and an elevation angle of about 
—61° and azimuthal angle of about 272° for the second one (FR2). The same chirality and 
similar axis orientations support the interpretation of a single ICME. Both flux ropes have 
strong axial magnetic held components compared with the azimuthal ones. Also note the 
largely southward orientation of the hux ropes. It is the strong axial held component on top 
of the largely southward hux-rope orientation, in conjunction with the relatively high solar 
wind speed, that may have resulted in the intense geomagnetic storm. 


4. Conclusions 

We have examined the sources of the 2015 March 17 and June 22 intense geomagnetic 
storms, the largest ones of solar cycle 24 up to the time of this writing. Key hndings 
are obtained on how the plasma and magnetic held characteristics of CMEs control the 
geomagnetic storm intensity and variability: 

1. A sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta scenario are proposed 
for the development of multi-step geomagnetic storms, based on the observed and recon¬ 
structed solar wind structures associated with the CMEs. The 2015 March 17 geomagnetic 
storm shows a two-step development, which is produced by the southward magnetic held 
components behind the preceding shock and those within two interacting CMEs. This falls 
into the sheath-ejecta-ejecta category. The 2015 June 22 geomagnetic storm exhibits a multi- 
step development, which is caused by the southward helds due to amplihcation by a series of 
preceding shocks and those within a single ejecta. This is classihed as a sheath-sheath-ejecta 
scenario. The multiple preceding shocks and sheaths may precondition the magnetosphere 
for the growth of an intense geomagnetic storm. 

2. We hnd two contrasting cases of how the CME flux-rope characteristics generate 
intense geomagnetic storms. Our GS reconstruction of the ejecta responsible for the 2015 
June 22 geomagnetic storm indicates that the geomagnetic storm resulted from the largely 





southward flux-rope orientation with a strong axial magnetic held component. However, for 
the 2015 March 17 geomagnetic storm the GS reconstruction reveals a much larger azimuthal 
held component than the axial component. The intense geomagnetic storm occurred despite 
low flux-rope inclinations. A southward flux-rope orientation is thus not a necessity for a 
strong geomagnetic storm to occur. 


3. The “perfect storm” scenario proposed bv Ihiu et al.l fj2014al) may not be as rare as 
the phrase implies. The 2015 March 17 intense geomagnetic storm occurred in spite of the 
relatively weak solar flares and an encounter with the CME flank. What makes it an intense 
geomagnetic storm is the interaction between two successive CMEs plus the compression by 
a high-speed stream from behind, which helps maintain strong ejecta magnetic fields and 
a relatively high speed. This is essentially the “perfect storm” scenario - a combination 
of circumstances results in an event of unusual magnitude, although the 2015 March 17 
event is not “super” in the same sense as the 2012 July 23 solar storm. Note that there are 
many combinations of circumstances that can occur to make an event more geo-effective, 
including pileup of events, pre-event rarefactions and held line stretching, shock enhancement 
of southward fields, and following high-speed streams causing compressions. This “perfect 
storm” scenario now seems useful and necessary to worry about because complex events with 
these combinations are common. 


The research was supported by the Recruitment Program of Global Experts of Ghina, 
NSFG under grant 41374173 and the Specialized Research Fund for State Key Laboratories 
of Ghina. We acknowledge the use of data from Wind and SOHO and the Dgt index from 
WDG in Kyoto. 
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Fig. 1.— Difference images of the 2015 March 15 (left) and June 21 (right) CMEs from 
LASCO C3 aboard SOHO. Note the interaction between the March 15 CME (CME2) and 
a preceding one that occurred on March 14 (CMEl). 
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March 2015 


Fig. 2.— Solar wind measurements at Wind and associated Dgt index for the 2015 March 
17 event. From top to bottom, the panels show the proton density, bulk speed, proton tem¬ 
perature, magnetic field strength and components, and D^t index, respectively. The shaded 
regions indicate two ICME intervals. The vertical dashed line marks the associated shock. 
The dotted curve in t he third paii el denotes the expected proton temperature calculated from 
the observed speed flLopea 119871) . The red and blue curves in the botto m p anel represent 


P^t va lues estimated using the formulae of lO’Brien fc McPherronI (120001) and iBurton et ah 
(119751) . respectively. 
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Fig. 3.— Reconstructed cross sections of ICMEl (upper) and ICME2 (lower). Black contours 
show the distribution of the vector potential, and the color shading indicates the value of 
the axial magnetic field. The location of the maximum axial field is indicated by the black 
dot. The dashed line marks the trajectory of the Wind spacecraft. The thin black arrows 
denote the direction and magnitude of the observed magnetic fields projected onto the cross 
section, and the thick colored arrows show the projected RTN directions. 
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Fig. 4.— Solar wind measurements at Wind and associated index for the 2015 June 22 
event. Similar to Figure 2. The shaded region shows the overall ejecta interval, while the 
hatched areas indicate two small flux ropes identihed within the ICME. Three shocks are 
observed ahead of the ejecta. The last shock (S4) was driven by the CME that occurred at 
the Sun on 2015 June 22 and was overtaking the ICME at 1 AU. 
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Fig. 5.— Reconstructed cross sections of FRl (upper) and FR2 (lower) identified inside the 
2015 June 22 ICME. Similar to Figure 3. These small flux ropes are both left-handed and 
have similar axis orientations. 






























